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Abstract. We explore the morphological and dynamical evolution of galaxy clusters in simu- 
lations using scalar and vector- valued Minkowski valuations and the concept of fundamental 
plane relations. In this context, three questions are of fundamental interest: 1. How does the 
average cluster morphology depend on the cosmological background model? 2. Is it possible 
to discriminate between different cosmological models using cluster substructure in a sta- 
tistically significant way? 3. How is the dynamical state of a cluster, especially its distance 
from a virial equilibrium, correlated to its visual substructure? To answer these questions, we 
quantify cluster substructure using a set of morphological order parameters constructed on 
the basis of the Minkowski valuations (MVs) . The dynamical state of a cluster is described 
using global cluster parameters: in certain spaces of such parameters fundamental band-like 
structures are forming indicating the emergence of a virial equilibrium. We find that the 
average distances from these fundamental structures are correlated to the average amount of 
cluster substructure for our cluster samples during the time evolution. Furthermore, signifi- 
cant differences show up between the high- and the low-f2 m models. We pay special attention 
to the redshift evolution of morphological characteristics and find large differences between 
the cosmological models even for higher redshifts. 

Key words. Galaxies: clusters: general - X-rays: galaxies: clusters - Methods: N-body sim- 
ulations - Methods: statistical. 



1. Introduction 

Galaxy clusters may be thought to constitute 
a sort of pocket guide to our Universe: al- 



background cosmology goes as follows: accord- 
ing to the hierarchical scenario, galaxy clus- 
ters were assembled through the merging of 
smaller obje cts, which collapsed first. Richstone] 



mological scales, they contain important mtor- 
mation about the Universe as a whole. One 
line of thought linking galaxy clusters and the 



tnougn tney are small m comparison to cos- et al ( 19 g^ suggested that the cluster dynam- 
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ical state is related to its age, which in turn de- 
pends on average on the present value of the cos- 
mological density parameter fl m . If, finally, the 
cluster dynamical state is mirrored by its sub- 
structure, one can establish a link between clus- 
ter morphology and the background cosmology 
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( "cosmology- morphology connection for galaxy 
clusters", Evrard et al. 1993). Therefore, the 



to a sample of Einstein IPC cluster images. 



cluster substructure may be a powerful tool to 
study the background cosmology. Summarising 
the results of the theoretical analyses (see also 



Bartelmann et al. 1993), one can state that in 



Buote & Tsai (1995) introduced the power ratio 
method, a technique based on the multipole ex- 
pansion of the two-dimensional potential gener- 
ating the observed surface X-ray brightness, see 
also|Buote fc Tsai (1996)); [Buote fc Xu (1997|); 



low-£l m cosmologies the clusters should on av- 
erage show a smaller amount of substructure 
than in high-f2 m models. Since this argument 
oversimplifies the complex dynamical situation 
in galaxy clusters, it has to be complemented 



Tsai fe Buote (1996|); |Valdarnini et al. (1999| ) 



Kolokotronis et al. (2001) studied the correla- 



tion between substructures observed both in the 
optical and X-ray bands. 

Cosmological N-body simulations have been 
used to test the dependence of cluster substruc- 



Note, that we need a thorough definition and ture on different cosmological models (Evrard 


description of cluster substructure for this ar- et al. 1993 




Mohr et al. 1995; Jing et al. 199£ ; 


gument 


L 


Thomas, P. A., Colberg, J. M., Couchman, H. 


In this context, it is still a difficult task to de- M. P., et al. 1998; 


West et al. 1988 


). Crone et al. 


scribe both the inner cluster state and the clus- (1996) applied different substructure statistics 



ter morphology quantitatively in a reliable way. 
- In this paper, therefore, we use new tools to 
quantify cluster substructure as well as the in- 
trinsic cluster state. We analyse cluster simula- 
tions with these tools and characterise the sub- 
structure of different cluster components, its re- 
lation to inner cluster properties and the differ- 



to galaxy clusters obtained in different cos- 
mological models from numerical simulations. 
They conclude that the "centre-of-mass shift" 
is a better indicator to distinguish between dif- 
ferent models than, e.g., the Dressier Shectman 
statistics (Dressier & Shectman 1988), which 



ences between cosmological background models 



as traced by the averaged cluster substructure. 
In particular, we test the theoretical assump- 
tions behind the "cosmology-morphology con- 
nection" . 

So far, various methods have been used to 
quantify the amount of substructure in galaxy 
clusters. In the optical band several tcch- 



does not provide significant results (Knebe & 
Mtillcr 2000| ). |Pinkncy et al. (1996p tested sev- 



niques ( 


Dressier & Shectman 1988 




Wcst & 


Bothun 1990|; 


Bird 1994 


) use the galaxy posi- 



tions and velocities. Other methods are based 



on the hierarchical clustering paradigm (Scrna 



fc Gerbal 1996| ; purzadyan fc Mazure 199g 
wavelet analysis (Girardi, M., Escalera, E 



Fadda, P., et al. 1997), or moments of the X-ray 
photon distribution ( Dutta 1995] ). 
X-ray images of galaxy clusters were also used 
to study substructure; contrary to optical clus- 
ters, they are scarcely contaminated by fore- 
and background effects. Mohr et al. (1995) ap- 
plied statistics based on the axial ratio and the 
centroid shift of isophotes ( Mohr et al. 1993| ) 



eral descriptors using N-body simulations and 
recommended a battery of morphology param- 
eters to balance the disadvantages of different 
methods. 

So far, however, a unifying framework for the 
morphological description of galaxy clusters is 
missing. Several aspects of cluster substructure 
have to be distinguished in order to provide 
an exhaustive characterisation. Also the connec- 
tion to a possible cluster equilibrium has not yet 
been scrutinised. 

In this paper we apply Minkowski valuations 
(MVs) ( |Mecke et al. 1994 |Beisbart et~al 
2001a| ,|b[) to cluster substructure and use fun 



From an observational point of view, there is 
clear evidence for the existence of substructure 



in galaxy clusters, both from optical data (Geller 


fc Beers 19821 


Dressier & Shectman 1988 


; West 


fc Bothun 1990 


Bird 1995) and from X-ray fm- 


ages (Jones & Forman 19921; Bohringer 1994 


; Mohr 


3t al. 1993 


)■ 



damental structures to quantify the dynamical 
state of galaxy clusters. The Minkowski frame- 
work provides mathematically solid and unify- 
ing morphometric concepts, which can be ap- 
plied to cluster data without any statistical 
presumptions. These measures distinguish ef- 
fectively between different aspects of substruc- 
ture and discriminate between different cos- 
mological background models. Our interest is 
both methodological and physical: on the one 
hand, we are looking for an appropriate method 
to quantify cluster substructure; on the other 
hand, we ask physical questions like: how are 
the Dark Matter (DM) and the gas distribution 
related to each other? 
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For our investigation, we employ combined N- 
body/hydrodynamic simulations. This simula- 
tion technique is particularly suitable for our 
purposes, since it traces both the dark matter 
and the gas component of a cluster. We con- 
struct relatively large data bases of cluster im- 
ages from the simulations which can be com- 
pared to real cluster images. 



spectral index of n = 1 and selected a baryon 
density parameter of fl^h 2 — 0.015. In the 
CHDM model we had one massive neutrino with 
mass m v — 4.65 eV, yielding a HDM density 
parameter £L n = 0.2. In the ACDM model the 
vacuum contribution to the energy density was 
f^A = 0.7 in accordance with recent observa- 



tions of Supcrnovae (Perlmutter, S., Aldering 



The plan of the paper is as follows: after an G., Goldhabcr, G., et al. 1999| ). Therefore, the 



explanation of the simulations and cosmologi- 
cal models in Sect. |[ we give an introduction 
into Minkowski valuations in Sect. ||. We em- 
ploy these tools in Sect. |J in order to compare 
the clusters within the different simulations. An 
analysis of fundamental plane relations is pre- 



dcnsity parameter of matter f2 m was 1 for CDM 
and CHDM, and 0.3 for ACDM. Since we are 
dealing with galaxy clusters and need a fair 
number of them, all models were normalised in 
order to match the present-day cluster abun- 
dance n c (M > M c ) = 4 ■ !Q- e h 3 for M c = 



sented in Sect. p[ We draw our conclusions in 4.2/i 1 10 14 M Q ( 


Eke et al. 1996; |Girardi, M., 


Sect. 


6. Escalera, E., Fadda, D., et al. 1997 


). Using these 



2. The cosmological models and the 
simulations 

In order to investigate cluster substructure in 
different cosmological models, a data base of 
galaxy clusters was generated on the base of 
TREESPH simulations. Three background cos- 
mologies were chosen differing both in terms of 
the values of the cosmological parameters and 
the power spectra. We restricted ourselves to 
CDM models; the simulations are described in 
more detail in Valdarnini et al. (1999| ), where 
also a morphological analysis was done using 
the power ratios (PRs, sec Buote fc Tsai 1995| ). 
We extend this work in several directions, e.g. 
by probing the morphological evolution and by 
connecting cluster substructure and inner dy- 
namical cluster state. 

Since observations indicate that the curvature 
parameter f^K vanishes (see for instance |de 



ernardis, P., Ade, P. A. R., Bock, J. J., et 



1. 200C ) we considered three spatially flat cos- 
mological models, namely two high-O m models 
(a standard Cold Dark Matter model - CDM - 
and a model where the Dark Matter consists of 
a mixture of massive neutrinos and Cold Dark 
Matter - CHDM) and one low-f2 m model (a 
model with a non-vanishing cosmological con- 
stant - ACDM). For the Hubble parameter we 
chose h = 0.5 for the CDM and the CHDM 
model, and h — 0.7 for the ACDM model; here, 
as usual, the Hubble constant is written in the 
form H Q = 100/ikms -1 Mpc -1 . 
With respect to the power spectra comprising 
the influence of the initial matter composition 
on structure formation we adopted a primeval 



normalisations only the ACDM model is consis- 
tent with the measured COBE quadrupole mo- 
ment at the la level. In order to reduce the influ- 
ence of cosmic variance the same random num- 
bers were used to set the initial conditions for 
all cosmological models. Therefore we roughly 
look at the same clusters in all cosmologies. 
The cluster simulation technique consisted of 
two steps: first for each model a large collision- 
less N-body simulation was performed using a 
P 3 M code in a box of length L = 200 a /i _1 Mpc, 
where a is the cosmological scale factor being 
one at present day. We considered N p = 10 6 par- 
ticles for the CDM and CHDM models, each, 
while N p = 84 3 particles were chosen for the 
ACDM model, the only low-density cosmology 
investigated here; thus the mass of one sim- 
ulation particle is approximately equal in all 
cosmological models. The simulations were run 
starting from an initial redshift z- lri , depending 
on the model (for more details see Valdarnini 
et al. 1999 ), down to z = 0. At the final redshift 



we identified galaxy clusters using a friend-of- 
fricnd algorithm in order to detect overdensities 
in excess of ~ 200 Sl^ 0,6 . For our further analy- 
ses, we took into account only the 40 most mas- 
sive clusters. 

As a second step we applied a multi-mass tech- 
nique ( |Katz fc White 1993j ; |Navarro et al. 1995| ): 
for each cluster we carried out a hydrodynamic 
TREESPH simulation in a smaller box start- 
ing from Zi n . For this we identified all cluster 
particles within r2oo (where the cluster den- 
sity is about 200 f^ ' 6 times the background 
density) at z — 0. These particles were back- 
tracked to Zi n in the original cosmological sim- 
ulation box. For each cluster a cube enclos- 
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Fig. 1. Two of the simulated cluster images to be analysed. They show the X-ray image of one 
cluster in the ACDM model at redshift zero (left panel) and the corresponding pseudo-DM-image, 
where the gas density has been replaced by the DM-density (right panel) . Both images are sampled 
within a circle of radius 1.5/i _1 Mpc around the peak of the gas/DM-surface brightness, respec- 
tively. Sec the main text for a further description of the image construction. 



ing all of these particles was constructed; its 
size L c was ranging from 15 to 25/i~ 1 Mpc. A 
higher-resolution lattice of Nt, — 22 3 grid points 
then was set into these cubes. Different lat- 
tices were used for the different mass compo- 
nents; to avoid singularities these lattices were 
shifted with respect to each other by 1/2 of 
the grid constant along each spatial direction. 
For the CHDM simulations the hot particles 



spatial resolution of the simulations can be es- 
timated by the ratio Egas/^ooj which never ex- 
ceeds a value of about 0.04. 
The numerical integrations were performed with 
a tolerance parameter 8 = 0.7 and using a leap- 
frog scheme for the time integration; the mini- 
mum time step allowed was 3 x 10 6 years for the 
gas particles and 6 x 10 6 years for the DM part. 



bear a small initial peculiar velocity following 



Katz 



a Fermi-Dirac distribution with Fermi velocity 
t>o = 5(1 + Zi„)(10eV/m y ) km s _1 . For the gas 
particles we started with an initial temperature 



Viscosity was treated as in Hernquist 
(1989| ) with a = 1 and = 2. The effects of 
heating and cooling were not considered in the 
simulations. Tests assessing the quality of the 



simulations are described in Valdarnini et al 



Tj = 10 4 K. The TREESPH simulation was then (1999| ). We saved numerical outputs at differ 



run using all particles which lie inside a sphere 
of radius L c around the centres of the cubes. 
The gravitational softening parameters e were 
the same for all clusters within each simulation 
and cosmological model. For the gas particles 
they were chosen to be e gas = 80, 100 , 60 kpc 
for the CDM, the CHDM , and most of the 
ACDM clusters, respectively. However, for the 
five most massive ACDM clusters e gas was set 
to 80 kpc. As softening parameters for the Dark 
Matter particles we took £d = 200, 231, 125 kpc 
for the CDM, CHDM, and ACDM model, re- 
spectively. For the simulation particles we ap- 

1/3 

plied the scaling e\ oc m ; . Note, that the soft- 
ening lengths were fixed within proper physi- 
cal space; however, the redshift Zi„ is chosen in 



ent redshifts, such that the cluster morphologi- 
cal evolutions could be investigated within the 
different models. 

Using the simulations we generated cluster im- 
ages which mimic observations in a realistic 
manner as follows: the gas density was esti- 
mated on a cubic grid with a grid constant 
of 0.03 /i _1 Mpc for each model. We took the 
square of this density at each grid point and 
calculated the approximate integral of p 2 along 
the line of sight orthogonal to a random plane 
(it is the same random plane for all clusters, 
simulations, and redshifts), with 101 x 101 pix- 
els. We considered the cluster as approximately 
isothermal, such that the X-ray emissivit y is just 
proportiona l to this integral (see, e.g., Tsai fc 



such a way, that the mean particle separation Buotc 1996[ ) 



is always smaller than the softening length. The We applied the same method also to the DM 



C. Beisbart et al.: Morphological evolution of galaxy clusters 



5 



particles; evidently, this does not lead to a phys- 
ically observable quantity. However, in this way 
we get the emissivity we would obtain if the gas 
distribution would trace the DM (a constant ra- 
tio between gas and DM distribution drops out 
in our analysis). We show both an X-ray and a 
DM-image in Fig. [y. The images are analysed 
using the Minkowski valuations, which are de- 
scribed in the next section. 



3. Minkowski valuations 

The Minkowski valuations (MVs) provide an el- 
egant and in a certain sense unique description 
of spatial data. They were introduced into cos- 



mology by Mecke et al. (1994) and have been 
applied to answer a number of questions re- 
garding the morphology of the large-scale struc- 
ture, see, e.g. Kerscher, M., Schmalzing, J. 



Retzlaff, J., et al. (1997|); 


Kerscher et al. (1998); 


Schmalzing & Gorski (199S ); 


Salmi ct al. (1998); 



Schmalzing, J., Buchert, T., Melott, A. L., et al 
(1999|); [Kerscher, M., Mecke, K. Schmalzing~ 



et al. (2001). So far, they were employed mainly 
in situations where perturbations of a homoge- 
neous background were to be expected and the 
amount of clustering had to be quantified. For 
galaxy clusters, however, the situation is differ- 
ent. Galaxy clusters are intrinsically inhomoge- 
neous systems, thus the main issue is how far 
their structure is away from a symmetric and 
substructure-poor state which does not show 
the influence of recent mergers. 
For this reason, additionally to the scalar 
Minkowski functionals, we use vector-valued 
Minkowski valuations (also known as "Quermafi 
vectors"), which feature directional informa- 
tionQ. In this section we give a short overview 
of both the scalar and the vector-valued MVs. 
For a general approach, let us consider pat- 
terns P, Q, i.e. compact sets within Euclidean 
space. A morphometric (geometrical and topo- 
logical) description of such spatial patterns is 
adequate, if it obeys a number of covariancc 
properties specifying how the descriptors change 
if the pattern is transformed. The Minkowski 
valuations are defined by three types of covari- 
ances. 



2 We reserve the name "Minkowski functionals" 
to the scalar Minkowski functionals. The entirety of 
both scalar and vector-valued Minkowski measures 
is referred to as "Minkowski valuations" . 



1. The first class of covariance refers to mo- 
tions in space: a morphological descriptor 
should behave in a well-defined way, if the 
pattern is moved around in space. The scalar 
Minkowski functionals are motion-invariant, 
whereas the vector-valued MVs transform 
like vectors (motion equivariance) . 

2. Secondly, the descriptors should obey a sim- 
ple rule specifying how they sum up for 
combined patterns whose set union is con- 
structed. Both classes of Minkowski valua- 
tions V obey the additivity: 



V(P U Q) = V(P) + V(Q) - V(P n Q) 



(1) 



3. Moreover, continuity states that a descriptor 
should change continuously if the pattern is 
distorted slightly^. 

These simple properties already define the MVs, 
since - at least on the convex rin g^ the Theorem 
of Hadwiger (Hadwiger 1957; [Hadwiger 



Schneider 1971 ) guarantees that in d dimen- 
sions, only (d + 1) of such measures - either 
scalar- or vector- valued ones - are linearly inde- 
pendent. 

Because of the structure of our cluster data, we 
concentrate on the case of d = 2. Here, the 
Minkowski functionals have intuitive meanings: 
Vo is the surface content of the pattern, V\ its 
length of circumference, and Vi its Euler char- 
acteristic \, which counts the components of 
the cluster and subtracts the number of holes. 
Note, that all of these functionals can be ex- 
pressed as an integral: Vq obviously is the two- 
dimensional volume integral of the pattern, V\ 
is its (one-dimensional) surface integral, and V% 
weights each surface element dS with the lo- 
cal curvature n. This integral representation is 
also valid for the vector-valued Minkowski val- 
uations; in this case, additionally, the integrals 
are weighted with the position vector; therefore, 
the Quermafi vectors are spatial moments of the 
scalar Minkowski functional^. Summarising we 



A closer analysis shows that this requirement 
can be only imposed for convex bodies. 

4 The convex ring comprises all finite unions of 
convex bodies. 

5 More precisely, they are first-order moments. 
Higher moments are invest igated as Minkowski t en- 
sors in forthcoming work (Beisbart et al. 2001b). 
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consider the following measures: 



V = / dV 



V 1 



AS 



V 1 



dVx 



(2) 



dP 



dSx 



dP 



V 2 = dSn , V 2 = dSnx , 

JdP JdP 

where k refers to the local curvature. - It is 
convenient to divide the vectors by the corre- 
sponding scalars to arrive at the (curvature) 
centroids^]: 



(3) 



The centroids localise individual morphological 
features. In particular, they coincide with the 
symmetry centre for spherically symmetric pat- 
terns. On the other hand, centroids constituting 
a finite triangle indicate the presence of asym- 
metry. 

Note, that the MVs obey covariance proper- 
ties with respect to a scaling of the pattern P, 
too: If we scale the pattern by a > to get 
ctP = {ax\x G P}, the scalar Minkowski func- 
tional transform like Vi(aP) = a d ~ l Vi(P); the 
vectors transform like Vj(<xP) = a d ~ t+1 V,(P). 
This is important since we sometimes have to 
compare data of different size. 
Obviously, cluster images are not patterns in 
the above sense, but rather consist of galaxy 
positions or pixelised maps reflecting the sur- 
face brightness within a certain energy band. 
Thus, we have to construct patterns from the 
cluster data. Here we use the excursion set ap- 
proach where we smooth the original data using 
a Gaussian kernel in order to construct a real- 
ization of a field u(x). The excursion sets 



M w = {x\u(x) > w} 



(4) 



are analysed with the aid of the Minkowski val- 
uations. Varying the density threshold w allows 
us to probe different regions of the cluster. The 
smoothing both reduces possible noise and picks 
out a scale of interest or resolution beneath 
which substructure is no longer resolved. 
The MVs, represented as functions of the den- 
sity threshold, contain very detailed informa- 
tion. In this paper, however, we want to com- 
pare cluster images drawn from a larger base 

6 In the following we speak of "centroids" or "cur- 
vature centroids", the latter term becomes more 
plausible in higher dimensions, where almost all cen- 
troids are connected to curvature. 



of simulated galaxy clusters. We are therefore 
interested in the average morphological cluster 
evolution in different cosmologies. In order to 
condense the detailed information present in the 
MVs, we construct robust structure functions 
which allow us to compare clusters of different 
size statistically. 

This can be done by integrating over the den- 
sity thresholds and weighting with functions of 
the Minkowski valuations. We define an average 
over different density thresholds via 



</>» 



1 



duf 



(5) 



Here, we consider three classes of robust struc- 
ture functions, which feature different aspects 
of the substructure and span a morphological 
phase space. We distinguish the following three 
classes of substructure, which are quantified us- 
ing one or more structure functions: 



1. 



the dumpiness C = J {(x — l) 2 ) u is a mea- 
sure of the number of subsystems in a clus- 
ter; 

the shape and asymmetry parame- 
ters {A = (Area(A(p 3 ))) u , A x = 
(perimeter(A(pj))) u ) refer to the de- 
gree of asymmetry and the global shape 
of the cluster ("is the cluster spherical or 
elongated?"); here we use the fact that 
curvature centroids which do not coincide 
within one point indicate the presence 
of asymmetry. In particular, the size of 
the triangle A(pj) which connects the 
curvature centroids (measured by its area 
and perimeter) serves as a measure of the 
asymmetry present within the cluster, 
the shift parameters Sj=o,i,2 = 
y/(\Pi - (Pt) u l 2 ) u account for the vari- 
ation or shift of morphological properties 
in a quantitative way by considering dif- 
ferent density thresholds (these parameters 
are generalisations of the frequently used 
centre-of-mass shift and centroid variation, 
see |Crone et al. (1996Q ; |Mohr et al. (1993| , 



An effective method to calculate Minkowski val- 
uations from pixelised maps can be developed in 



7 It turns out that S2 is closely related to the 
dumpiness. - For the sequel we concentrate on the 
structure functions C, Ai, and Si. We successfully 
tested the robustness of these latter structure func- 
tions. 
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analogy to Bchmalzing & Buchcrt (1997) and is 
given in 



Beisbart, C, et al. (2001) 



4. Cluster substructure and the 
background cosmology 

To start our analysis of the simulated clus- 
ters, we probe the connection between the back- 
ground cosmology and the cluster morphology. 
Since in this case not so much the substructure 
of individual clusters rather than the mean mor- 
phology is of interest, we define cluster samples 
consisting of all clusters at one redshift within 
one model (unless otherwise stated we analyse 
one random projection per cluster)^]. In order to 
trace the mean substructure evolution, we av- 
erage the structure functions over all clusters 
within one sample. 

What is the amount of cluster substructure, and 
how does it evolve within the three cosmologi- 
cal models? - The simulated clusters are "ob- 
served" within a quadratic window of 3/i _1 Mpc 
width centred at the peak position of the sur- 
face brightness. The data are smoothed us- 
ing a Gaussian kernel with different smoothing 
lengths in order to reduce the sensitivity to noise 
and to probe different scales of the substruc- 
ture. We concentrate on intermediate values of 
the smoothing scale A (A ~ 0.05 — 0.15/i _1 Mpc, 
Cen ( 1997| ) employs values of the same order of 
magnitude )[]. To define the cluster on the im- 
age, we draw circles around the peak with radii 
r w = 0.8, 1.0, 1.2, l^/i^Mpc (this definition is 
in the spirit of Abell's cluster identification in 
the optical, see Abell 1958| , we call this circle 
the cluster window) and neglect the rest. The 
integration limit w m i n in Eq. (|^) is chosen to 
be twelve times the background which is deter- 
mi ned from the rest of the image similarly as 



in |Bohringcr, H., Vogcs, W., Huchra, J., et al 



(2000 ) , u max in Eq. (|5|) is the maximum cluster 
surface brightness. 

In Fig. H we show results for the X-ray clus- 
ter morphology within the three models hav- 
ing applied a smoothing length of 0.05 /i _1 Mpc. 

8 On account of applications below, we discarded 
a few clusters which showed obvious pathologies 
such as a strong bimodality. We have 39 clusters 
for the CDM model, 37 for the CHDM and 35 for 
the ACDM model. 

9 A smoothing length of 0.05fr _1 Mpc is smaller 
or equal than the gravitational softening length for 
the gas and not much smaller than that one for the 
DM. 



A couple of things are obvious at first glance: 
there is a significant difference between the high- 
ly models (CDM and CHDM) and the low- 
Sl m model (ACDM). These differences are vis- 
ible in most of the structure functions and are 
in accordance with the theoretical expectations: 
the low-f2 m model shows by far less substruc- 
ture than the other two models - at least for 
most redshifts investigated here. The CDM and 
CHDM models, however, do not seem to be 
distinguished well. Therefore, the morphology- 
cosmology connection is mainly sensitive to the 
values of the cosmological parameters, but per- 
forms poorly in discriminating between differ- 
ent power spectra. The dumpiness is particu- 
larly sensitive. - Regarding the redshift evolu- 
tion, a clear trend is visible towards more re- 
laxed and substructure- poor clusters. In partic- 
ular, there are also large morphological differ- 
ences between the models at higher redshifts. 
Morphological evolution of galaxy clusters may 
therefore serve as a more sensitive test than the 
present-day cluster morphology. Note, that the 
averaged morphology evolution still looks rela- 
tively spiky. The reason is that for individual 
clusters the evolution of the structure functions 
proceeds in a discontinuous manner, when sub- 
clumps enter the cluster window. Therefore, one 
has to average over several clusters in order to 
get a typical morphological cluster evolution. 
The structure functions A.; and Si have a dimen- 
sion and therefore quantify the absolute amount 
of substructure. To investigate the substructure 
relative to cluster size, we normalise Ai and 
Si to the individual cluster size estimated via 
the two-dimensional half-light radius around the 
peak of the X-ray surface brightness. As visible 
from the bottom row of Fig. ||, the qualitative 
evolution and the differences between the mod- 
els are similar as before. 

So far we concentrated on the morphological 
evolution as traced by the X-ray luminosity and 
thus the X-ray gas. But is also the DM morphol- 
ogy different for the cosmological models? The 
results in Fig. || show the mean substructure 
evolution for galaxy clusters (A = 0.05/i -1 Mpc) 
and indicate that the DM morphology in clus- 
ters is even more sensitive to the cosmological 
background than the gas. 

In order to strengthen our claims and to com- 
pare the performance of the gas- and the DM- 
morphology in a systematic way, we take into 
account the whole distribution of the struc- 
ture functions for our cluster samples. The 
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Fig. 2. The averaged morphological evolution of the galaxy clusters in our three cosmological 
models. The ensemble-averaged structure functions dumpiness, shape and asymmetry, and shift 
of morphological properties were determined from the X-ray luminosity for a smoothing length of 
A = 0.05/j _1 Mpc and are shown vs. rcdshift (first row). CDM: solid line; CHDM: short dashed 
line; ACDM: long dashed line. We consider a spherical window with radius r w = 1.4/i -1 Mpc; C is 
dimcnsionless; Ai and Si are given in units of (/i _1 Mpc) 2 and ft. _1 Mpc, respectively. In the second 
row we show the A^ and S^, where Ai and Si were scaled to cluster size (an estimate of the half 
light radius) in order to get dimensionless parameters. 
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Fig. 3. The morphological evolution of the DM within the clusters. The ensemble-averaged 
structure functions dumpiness, shape and asymmetry, and shift of morphological properties for 
a smoothing length of 0.05ft. _1 Mpc are shown vs. redshift. CDM: solid line; CHDM: short dashed 
line; ACDM: long dashed line. Contrary to Fig. here we probe the DM. We consider a spherical 
window of radius r w = 1.4/i _1 Mpc; for the units see the caption of Fig. ^. This time, we also 
show the one-sigma fluctuations around the average structure functions for the ACDM model. 
Note, that large fluctuations are to be expected on physical grounds because of the wide range of 
clusters investigated, having different formation times and environments. We conclude that one 
has to consider the whole morphology distribution within the samples in order to get significant 
claims regarding the background cosmology (see below). 
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Kolmogorov-Smirnow test is a suitable tool to 
answer the question whether two data samples 
are likely to be drawn from the same distribu- 
tion. It measures the distance between two cu- 
mulative distributions D\(X) = p(x < X) and 
Z?2 (estimated from the empirical data) via 



d KS = max.{\D 1 (x) 



D 2 (x)\} 



(6) 



For any value of c?ks one can estimate the prob- 
ability pks that the same random process gen- 
erates two samples being "farther away" from 
each other than (Iks- If ^ks is large and, ac- 
cordingly, pks is small, the distributions under 
investigation are likely to stem from two differ- 
ent ensembles. Table [j] presents results of a KS 
test for redshift z = and a smoothing length of 
A = 0.05fo Mpc. The small values of j?ks show 
that the distinction of the cosmological models 
by means of the structure functions is effective 
regarding the values of the cosmological param- 
eters. A systematic imprint of the power spec- 
trum on cluster morphology, however, does not 
seem to be noticeable. The DM is even better 
than the gas in discriminating between the back- 
ground cosmologies. For instance, in comparing 
the CDM and the ACDM model, the probabil- 
ity of the null hypothesis is smaller than 10 -5 
using the dumpiness (A = 0.1/i _1 Mpc). 
So far, we investigated only one smoothing scale 
and one size of the cluster window; but one may 
ask which cluster regions are most interesting 
for a distinction between cosmological models 
and which smoothing lengths are optimal for 
our purposes. In order to answer these ques- 
tions, we first focus on the gas and estimate our 
structure functions for all cluster samples at red- 
shift z = for a number of smoothing lengths 
and scales of the cluster window. For each win- 
dow and smoothing scale we calculate the KS- 
distances between our three models. The results 
show that for the dumpiness small smoothing 
lengths A are more favourable than larger ones. 
For lower resolutions, therefore, the subclumps 
are smeared out, and the dumpiness is domi- 
nated by random fluctuations. 
On the other hand, the discriminative power 
of the dumpiness is enhanced for larger scales 
of the window. The reason is that subclumps 
which have not yet merged with the main clus- 
ter component are to be found at the outer clus- 
ter parts. The other structure functions mostly 
depend only relatively weakly on the window 
scale. We conclude that the outer cluster re- 



gions, which probably have not yet been viri- 
alised, are of more interest for the cosmologist. 
Moreover, the values of the structure functions 
rise for larger windows. This confirms results 
by Valdarnini et al. (1999), who found a similar 
behaviour (see, e.g., their Table 5). For the DM 
morphology one obtains comparable results. 

5. Fundamental plane relations 

Using cluster simulations one can test the basic 
assumption behind the morphology-cosmology 
connection presuming that the morphology 
of a cluster mirrors its inner dynamical state 
reliably. In this section we try to bring together 
morphology and inner state of our galaxy 
clusters. Mostly, we focus on the ACDM model 
as the nowadays favoured one. 
Observationally there is evidence that clusters 
of galaxies undergo a dynamical evolution lead- 
ing to a sort of equilibrium. This equilibrium 
seems to be manifest in fundamental plane 
relations holding within three-dimensional 
spaces of global cluster parameters where 
clusters tend to populate a plane. Since the 
cluster parameters are logarithms of observable 
quantities, the fundamental plane (FP) corre- 
sponds to a power law constraint among the 
real cluster parameters. Usually fundamental 
plane relations are explained in terms of the 



virial theorem of ( Chandrasekhar & Lee 1968) 



which, however, is strictly val id only for isolated 
syste ms (for a discussion see Fritsch fc Buchert 
19991) . 



There are several interesting spaces of global 
cluster parameters, depending on whether 
optical or X-ray data are available. Usually, 
the scale of the cluster, an estimate of its mass 
and a quantity related to its kinetic energy 
like the velocity dispersion of the galaxies or 
the temperature of the X-ray emitting gas are 
considered, see, e.g., Bchaeffcr ct al. (1993[ ); 



Adami ct al. (1998) for optical fundamental 
planes and Annis (1994 ); [Fritsch fe Buchert] 
(I999| ); [Fujita fc Takahara (1999j ) for X-ray 



clusters. In each case, indirectly, the potential 
and the kinetic energy are referred to. 
Fritsch fc Buchert (1999[ ) showed that the 
substructure of a cluster is correlated to its 
distance from the fundamental plane using the 
COSMOS/ APM and the ROSAT data. In the 
spirit of their work, we try to establish a similar 
connection for simulated X-ray clusters. 
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Table 1. The discriminative power of the structure functions regarding the cosmological back- 
ground models. The X-ray morphology is considered for a smoothing length of A = 0.05/i~ 1 Mpc 
and at redshift zero. For each pair of background models, we calculate the Kolmogorov-Smirnow 
distance c?ks between the distributions of the structure functions as well as the probability of the 
nullhypothesis. As one can see, the gas substructure clearly discriminates between the high- and 
the low-17 m models. 





Gas, 


r w = 1.40ft" 


1 Mpc, z 


= 0, A = 


0.05ft _1 Mpc 




CDM 


- ACDM 


CHDM 


- ACDM 


CDM - CHDM 






PKS 


dKS 


PKS 


dKS PKS 


c 


0.29 


7.7 % 


0.31 


4.5 % 


0.17 57.7 % 


Ax 


0.35 


1.6 % 


0.40 


0.4 % 


0.12 91.4 % 
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5.1. Fundamental plane relations for the 
simulated clusters 

Using our simulated clusters, we test three pos- 
sible parameter spaces spanned by: 1. the mass, 
the half-light-radius, and the emission- weighted 
X-ray temperature; 2. the mass, the half-mass- 
radius, and the X-ray luminosity; 3. the mass, 
the half-mass-radius, and the velocity disper- 
sion. In all of those global parameter spaces, one 
can observe a band-like fundamental structure. 
This thin band may be fitted either by a plane 
or a line. In Table || the global cluster param- 
eters are summarised using a compact notation 
which we will use from now on. 
The parameters defining the different cluster 
parameter spaces are estimated from the sim- 
ulations as follows: the cluster mass is quan- 
tified via M200 contained within an overden- 
sity 8 C times the critical density p c : M 2 oo = 
if ^PcHjoo) where S c ~ 178 ■ O" 45 in a flat 
cosmology ( poles fc Lucchin 1994 ) and where 
r2oo is the size of this overdensity. r2oo as well 
as rh are determined from the three-dimensional 
mass distribution around the density maximum, 
rh is the half-mass radius. T cm denotes the 
emission-weighted temperature of the gas, cal- 
culated from the gas thermal energy assuming 
an ideal gas. The bolometric X-ray luminosity 
is defined as L x = J dV(j^) 2 A c , where p gas 
is the gas density, /i = 0.6 the mean molecular 
weight, mp the proton mass, and A c the cooling 
function. In order to perform the volume inte- 
gration, the standard SPH estimator has been 
applied (Navarro et al. 1995), the summation in- 
cludes all particles within the virial radius f2oo- 
The velocity dispersion is estimated from all 
types of simulation particles. 
Before investigating the relationships between 



these parameters, we ask whether the distribu- 
tions of these parameters are consistent with 
each other for the different cosmological mod- 
els (in the sense of the KS test). We find consis- 
tency apart from the luminosity (which is higher 
on average for the ACDM model) and the half- 
light radius (clusters seem to be more compact 
within the CHDM model). 
We fit a plane and a line to each cluster sam- 
ple separately (one sample means one cosmo- 
logical model at one redshif t) using an ortho- 
gona l distance regression, see Boggs et al. (1987 , 
1989] ); this technique treats all variables the 
same way, i.e. no parameter is a priori thought 
of as dependent on the others^ The planes are 
parametrised by 



py 



ith plane 



(FP i) 



the lines are defined via: 



ith line 



(FL i) 



We call the best-fitting planes and lines funda- 
mental planes/ fundamental lines, respectively. 
We show one of the fundamental structures in 
Fig. I at redshift z = for the ACDM model. 
To get a clearer representation, we fit a sec- 
ond plane to our data under the constraint, 
that it be orthogonal to the fundamental plane, 
Takahara (1999|) did. A second 



as 



Fujita 



10 Note, that we do not lump together different 
models or redshifts for the fittings. In part, we used 
the ODRPACK software package for the orthogonal 
distance regression. 
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Table 2. Summary of the cluster parameter spaces investigated: the ith parameter space is 
spanned by the parameters P/ — ' ' for i 



1,..,3. 
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Fig. 4. The fundamental structure in the second parameter space as defined in Table 0. We 
consider the ACDM model at redshift zero. The first row shows the dependencies among the 
logarithmised physical parameters for all clusters. In the second row, we give a visual impression 
of the fundamental structure. The coordinate axes x\, X2 and £3 are determined in such a way, 
that the FP coincides with the x\ — x^ plane. 



constraint equation among the cluster param- 
eters should force the clusters to lie on a line 
in the global parameter space, this line should 
lie (more or less) on the intersection of both 
planes. We define a rotated coordinate system 
{xi, X2, £3} in such a way that the first (i.e. 
fundamental) plane coincides with the x\ — x-i 
plane, and the best-fitting orthogonal plane lies 
within the X2 — X3 plane. In this coordinate sys- 
tem, the scatters around both planes are eas- 
ily discernible as the 23- and cci-values for the 
clusters. The morphology of the structure obvi- 
ously is more band- than plane-like confirming 
results by Fujita fc Takahara (1998| ) , who call 



the structure they find in a different parame- 
ter space "the fundamental band" . This is also 
true for the other parameter spaces. A visual 
inspection of the fundamental structures shows 
furthermore that most outliers, which tend to 
prolongate the line, wander towards the bulk 
for lower redshifts. Note, that in our analysis 
statistical outliers are not removed. 
The values of the best-fit parameters are listed 
in Table || together with their 95%-confidence 
intervals. Since we have no measurement errors 
for the global parameters, we can give error bars 
only assuming the goodness of the fit. In order 
to probe whether our fundamental planes may 
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of the cluster and / denotes the baryon frac- 
tion. Assuming furthermore that Ri oc and 
M 2 oo oc M, we derive from our first fundamental 
plane-fit (for the ACDM model at redshift zero; 
we assume that the additional scaling relations 
used do not introduce additional uncertainties) 



fM „ T .1.61±0.24 T ,1.29±0.15 

jM 200 <xr h 1 



(7) 



Table 3. The best-fitting fundamental planes 
together with the 95%-confidence regions for all 
models and all parameter spaces considered at 
redshift zero. 



This result is consistent with theirs provided 
that the baryon fraction depends relatively 
strongly on the temperature: / oc T 73 . For 
comparison, using the same observational data 



as [Fujita fc Takahara (1999| ), [Mohr et al. (1999|) 

fmd //cm oc T u - 34±u ^ within r 500 , see |Mohi 



explain observed fundamental plane relations, 
we compare the parameters with simple the- 
oretical scaling laws and observed parameters. 
The problem about such comparisons, however, 
is that the definitions of the cluster parameters 
significantly depend on the techniques used to 
determine them from observations. Therefore, 
we have to use additional assumptions; we con- 
strain ourselves to X-ray fundamental planes 
at redshift zero. - From a theoretical point of 
view the third fundamental plane is the sim- 
plest one. A virial equilibrium requires that 
M V i r oc Ra 1 for the whole mass M, the scale 
R and the velocity dispersion cr 2 of a cluster 
being in virial equilibrium. If we simply iden- 
tify these parameters with the quantities span- 
ning the third cluster parameter space, we see 
that our values for /?| are consistent with the 
virial equilibrium for all models; moreover, /3| 
is compatible with the virial prediction for the 
ACDM model, and marginally consistent within 
the CDM model, but inconsistent for the CHDM 
model. A physical reason may be that, because 
of the high value of S7 m , a virial equilibrium 
is not yet reached for most clusters within the 
CHDM model. Perhaps also the plane-fit is de- 
termined by a few clusters not yet in equilib- 
rium; but certainly larger cluster samples are 
required in order to clarify this point definitely. 
The first fundamental plane can be compared to 
the results by Fujita fe Takahara (1999| ). Using 
data from Mohr et al. (1999) they find that the 



central gas density p g $ oc R± ' T enl , where 
Ri is the core radius. In order to relate our pa- 
rameters to theirs we estimate p 9i o by p g .o oc 
Mg as /r^ oc fM/r^, where M is the whole mass 



et al. (1999) for the exact definitions of the 
quantities they uscPj For the other cosmolog- 
ical models, the agreement is better. - From 
a theoretical point of view, one would expect 
that for a hydrostatic and a virial equilibrium 
M oc RT. For our data the T-dependence is 
slightly weaker, whereas R has a stronger influ- 
ence on the mass. 

For the second param eter space, we can use re- 
su lts by Fritsch (1996 ) whic h constitute the base 
of Fritsch fc Buchert (1999| ) . Assuming that the 
mass-to-light ratio is constant for galaxy clus- 
ters without scatter, our second fundamental 
plane translates into 



T ~ ^1.60±0.17 r 0.29±0.05 

l, oc r h 



(8) 



where L Q is the optical luminosity. On the other 
hand, putting together the virial mass estimat e 
and the fundamental plane from Fritsch (1996 ), 
relating L , Lx, and the optical half- light-radius 
r Q , we get 



=oc Li 



(9) 



Especially the dependence on the scale is consid- 
erably stronger in our fundamental plane, but 
again the discrepancies may be explained by the 
fact that our estimates of the cluster mass and 
scale differ from Fritsch' ones. 
Altogether, our results are in rough agreement 
with most of the theoretical expectations and 
the observed scaling laws. In detail, however, 
there are some inconsistencies to be found; but 
these incompabilities may be explained either 
with statistical fluctuations or by questioning 
some of the assumptions used in order to relate 
parameters estimated in different ways. 



11 Note, that [Vlohr et al. (1999[ ) use 90 % confi- 
dence regions instead our 95 % confidence levels. 
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Fig. 5. For these plots planes were fitted to the fundamental structures in each parameter space. 
We show the mean scatters around these fundamental planes (lower curves) and around the best- 
fitting orthogonal planes (upper curves) for each parameter space as a function of redshift. Again, 
we have: CDM: solid line; CHDM: short dashed line; ACDM: long dashed line. For technical 
reasons, we consider the CHDM model starting from z ~ 0.33, only. One sees that the scatter 
around the fundamental planes is clearly smaller than for the orthogonal planes. 
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Fig. 6. The exponents determining the second fundamental plane with their errors at low redshifts, 
see Eq. (FP i). Since we do not have measurement-like errors, the 95% confidence regions visible 
in the plot are estimated assuming the goodness of all fits. CDM: solid line; CHDM: short dashed 
line; ACDM: long dashed line. 



To analyse the morphologies of the funda- 
mental structures and their redshift evolutions 
quantitatively, we investigate the mean scat- 
ters around the fundamental planes, ctfp = 

\J Jq SiLi , and the orthogonal planes, ctop = 

\J Sili ^ii where we sum up the quadratic 
distances of the N clusters from the fundamen- 
tal planes, di, and the orthogonal planes, d{. As 
one can see from Fig. |[ erpp is decreasing on 
the whole for each of the first two global param- 
eter spaces. This, however, is not valid for the 
third parameter space if one fits the fundamen- 
tal structure using a plane. These details may 
indicate, that the band in the third parameter 
space is better fitted using a line. This conclu- 
sion is confirmed if one takes into account the 



scatter around the orthogonal plane, cop: Fig- §| 
shows that for the third parameter space, the 
scatter around the fundamental plane is only 
about two times larger than that one around 
the orthogonal plane. 

The evolution of one set of FP parameters is 
shown in Fig. || for low redshifts. The cosmo- 
logical models' confidence intervals, which were 
estimated again assuming the goodness of the 
fit, overlap for small redshifts (z < 0.05) indi- 
cating the consistency of the models regarding 
the location of the second fundamental plane. 
Apart from the CDM model the FP-exponents 
do not show any significant evolution for red- 
shifts z < 0.15. Similar results hold for the first 
parameter space. 

The scatters around the best-fitting lines are de- 
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Fig. 7. Now the fundamental structures are modelled using a line. The mean scatters around the 
best-fitting lines are shown for each model as function of redshift. Linestyles as in Fig. ||. 



creasing as a function of redshift in most cases 
(Fig. ^). Especially, the first and the third pa- 
rameter space show a strong redshift evolution, 
whereas for the second parameter space results 
are less definitive. This complements our ear- 
lier observations, that within the second param- 
eter space the fundamental structure is more 
plane-like, whereas the third fundamental struc- 
ture resembles a narrow band. To summarise the 
properties of the fundamental structures: in the 
first global parameter space we see a band-like 
structure, the structure in the second space can 
be understood as a plane, whereas in the third 
space the data are better fitted to a line. Using 
the corresponding fittings, the scatters go down, 
which we interpret in terms of an equilibrium 
attracting the clusters. 

5.2. Fundamental structures and cluster 
morphology 

So far our results indicate that the galaxy clus- 
ters are attracted by a quasi- equilibrium state 
mirrored by fundamental structures which seem 
to be more or less universal for all cosmolog- 
ical models. The evolutions of the mean scat- 
ters show that the clusters are approaching this 
quasi-equilibrium state in time in a sort of re- 
laxation process. In order to quantify how far in- 
dividual clusters are away from this equilibrium 
state, one can estimate their distances from the 
fundamental plane within each of the parame- 
ter spaces. The concept of a distance within the 
parameter space thus allows us to measure the 
inner dynamical state of a cluster. 
The physical nature of this quasi-equilibrium 
state can be confirmed if one can show that dif- 
ferent global cluster characteristics are accom- 



panying this evolution. An excellent candidate 
is cluster substructure; indeed, the cosmology- 
morphology connection is based on the assump- 
tion that the age of a cluster and therefore its 
dynamical state is reflected by cluster substruc- 
ture. We have already seen that on average, 
both the cluster substructure and the distances 
from the fundamental planes are decreasing in 
time. 

Therefore, we ask in a first step whether the 
sample-averaged substructure and the sample- 
averaged scatter around the fundamental struc- 
tures are correlated during their time-evolution. 
A basic test relies on Kendall's r, a non- 
paramctrical correlation 
193S| : ~ 



coefficient (Kendall 
Fritsch & Buchert 1999). In general, the 



amount of r reflects the strength of the corre- 
lations between two quantities within a given 
data set, while the sign of r specifies whether 
positive or negative correlations hold among the 
data points. Only values of r where p(r) (the 
probability of the nullhypothesis that no corre- 
lations among the data points exist) is smaller 
than 0.05 evince a statistically significant cor- 
relation. The results shown in Table || indicate 
that, for the case of the line- fitting, strong cor- 
relations exist for the CDM and for the ACDM 
model. For the CHDM model we have fewer red- 
shifts, a fact, that in part may explain the less 
meaningful results. 

To become more specific, we ask in a second 
step whether the distance from the fundamen- 
tal plane and the cluster morphology, quantified 
by the structure functions, are connected for 
individual clusters. Tentatively we carried out 
Kendall tests for each cluster sample relating 
the substructure parameters and the distances 
from the fundamental structures. Since we con- 
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Table 4. Correlations (Kendall's r) between the averaged substructure (as measured by the 
structure functions C, A 1; and Si for r w = 1.4 and a smoothing length of A = 0.05/i _1 Mpc) 
and the mean quadratic scatter around the fundamental lines. We consider all models and all 
parameter spaces. Similar results can be obtained using the plane-fitting, whenever the scatter 
around the plane is decreasing with time. The fact, that almost no significant correlations occur 
for the CHDM model, may be explained in part by noticing that we use fewer pairs of data points 
for the correlation analysis in this case. 
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Fig. 8. We apply a Kendall's test to relate the substructure and the dynamics of individual clusters. 
The substructure is measured using the dumpiness (estimated after having employed a smoothing 
length of 0.05/i~ 1 Mpc), the dynamical state of the cluster is quantified with the distance from the 
first fundamental line. We show the results, the Kendall correlation coefficient, as a function of 
redshift for all models. Note, that for significant results, we need a r larger than about 0.2 (slightly 
depending on the number of clusters which is different for the different models.) 



sidered several structure functions at different 
values of the smoothing scale, there is quite a 
lot of freedom. However, the results, as shown, 
for example, in Fig. |j, do not establish a signif- 
icant correlation between the substructure and 
the distance from the fundamental plane for in- 
dividual clusters; there is no connection persist- 
ing in time and throughout all of the cosmo- 
logical models. More specifically, there are not 
many significant correlations to be found at all; 



and some of them even turn out to be anticor- 
relations meaning that substructure-poor clus- 
ters are farer away from the fundamental plane 
than the substructure-rich ones. But there are 
also a number of positive correlations between 
substructure and distance from an equilibrium 
to be discovered for other models at certain red- 
shifts. 

The lack of a statistical significant relation be- 
tween substructure and dynamics for individual 
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clusters may have several reasons; in particular, 
a physical connection may be obscured on sta- 
tistical grounds. For example, if the fit of the 
fundamental plane is determined by a few clus- 
ters far from equilibrium, then the fitted funda- 
mental structure is distorted with respect to the 
real one and the distances from the fundamen- 
tal plane become distorted as well. This can be 
clarified in future analyses with larger cluster 
samples. 

It is, however, worth noticing that anticorrela- 
tions frequently occur in cases where the mean 
scatter around our fits of the fundamental struc- 
ture fails to decrease during the dynamical evo- 
lution. In other words, careful and appropriate 
fittings which result in a decreasing scatter with 
time reduce in part the indefinitive results in 
favour of a positive connection between sub- 
structure and the distance from the fundamen- 
tal plane. 

The results on the fundamental structures can 
be summarised as follows: 

1 . There is clear evidence for the existence of a 
fundamental band-like structure in all sorts 
of parameter spaces investigated here, even 
for moderate redshifts (< 0.4). This band 
can be fitted either using a plane or a line. 

2. Fitting the fundamental structures appro- 
priately, one can observe that in most cases 
the scatter around these structures decreases 
with time on the whole. 

3. Examining the evolution of the scatter 
around the fundamental structures and their 
morphologies one can specify, that in the 
first and third parameter space a line can 
be seen, whereas the data within the second 
space are better fitted by a plane. 

4. There are no large differences at redshift zero 
among the fundamental planes observed in 
different cosmological models. 

5. The scatters around the fundamental struc- 
tures are comparable for all cosmological 
background models investigated here0. 

6. There are some positive correlations between 
the structure functions and the distances 

12 This can be seen from Figs. |H| and |^. A quan- 
titative analysis based on the KS test shows that 
indeed the cosmological models are in most cases 
compatible with each other regarding the distances 
from the fundamental lines/planes for low redshifts. 
Thus, it is not possible to discriminate between the 
different background models merely by means of the 
scatter around the fundamental structures. 



from the fundamental structures (if appro- 
priately fitted) for the averaged cluster evo- 
lution; for individual clusters, however, the 
results are not conclusive. 

6. Conclusions 

The shape of a body frequently is its first prop- 
erty to be recognised. However, when the body 
is investigated more closely, more than a quali- 
tative description of its morphology is required. 
One wants to know, to what extent its appear- 
ance reflects inner properties, and one wants 
to compare the shape of the body to predic- 
tions of analytical models. For both purposes 
one needs robust descriptors of the body's mor- 
phology. This is also true for galaxy clusters. In 
this paper we investigated both the morphology 
and the inner dynamical state of galaxy clusters 
using large samples of simulated galaxy clusters. 
In order to measure the cluster substructure 
we employed structure functions based on the 
Minkowski valuations; the inner cluster state 
was quantified by the distance from a fundamen- 
tal band-like structure observed in a parameter 
space of global cluster descriptors. The inten- 
tion of our paper was twofold: we first tested 
a new method to measure the substructure of 
galaxy clusters. On the other hand, we investi- 
gated how far the cluster inner dynamical state 
is mirrored by the morphology and how differ- 
ent cosmological background models are distin- 
guished by means of the cluster substructure. 
Regarding the morphometry of galaxy clus- 
ters, i.e. the quantitative description of 
their size, shape, connectivity, and symmetry, 
the Minkowski functionals together with the 
Quermafi vectors allow for a discriminative and 
complete characterisation. They are based on a 
number of covariance properties and thus rest 
on a solid mathematical basis. The structure 
functions constructed from the MVs feature dif- 
ferent aspects of substructure successfully. 
Employing these methods we showed that the 
substructure of X-ray clusters distinguishes be- 
tween cosmological models in an effective way. 
As expected theoretically, the substructure on 
the whole is minimal for the ACDM model and 
higher for the high-f2 m models considered here. 
The power spectrum does not seem to have a 
systematic influence. We mainly focused on sim- 
ulated X-ray images; but also the DM substruc- 
ture can distinguish between the cosmological 
models. 
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Another important issue is the connection be- 
tween substructure and fundamental plane re- 
lations. This connection has not been investi- 
gated so far using numerical N-body simula- 
tions. In general, the evolution of fundamental 
plane relations within N-body simulations has 
not yet been scrutinised extensively. We could 
show that there are stable fundamental band- 
like structures within most cosmological models. 
Moreover, we found a positive correlation be- 
tween the averaged distance from this structure 
(if it is fitted appropriately) and the sample- 
averaged structure functions during time for two 
of our cosmological models. For individual clus- 
ters, however, we failed to produce definitive re- 
sults. Further investigations using larger simula- 
tions are in order to tackle this point. However, 
altogether there are weak indications that both 
our structure functions feature those aspects of 
substructure that reflect the inner cluster state 
and that the fundamental structures are the im- 
print of a physical equilibrium. 
These results raise a couple of new questions: 
how can we explain the fundamental structures? 
What is the physical origin of the degener- 
ated fundamental line? Are the fundamental 
bands dependent on the environment as sug- 
gested by [Miller et al. (1999| )? What is the pre- 
cise time evolution of fundamental structures? 
A number of tasks still remain to be done: in 
this paper, the FP-parameters were defined us- 
ing the three-dimensional clusters. How signif- 
icant are all the results found here, when one 
moves to more observation-like defined quanti- 
ties? In our results one thing is paradoxical: on 
the one hand, the mean substructure discrimi- 
nates well between the models, whereas on the 
other hand the mean scatters around the fun- 
damental bands are comparable for all sort of 
models. We conclude that the morphology is re- 
ally necessary to establish a connection between 
clusters and the global cosmological parameters. 
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